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Abstract. We introduce a <9-formulation of the orthogonal polynomi- 
als on the complex plane, and hence of the related normal matrix model, 
which is expected to play the same role as the Riemann-Hilbert for- 
malism in the theory of orthogonal polynomials on the line and for the 
related Hermitian model. We propose an analog of Deift-Kriecherbauer- 
McLaughlin-Venakides-Zhou asymptotic method for the analysis of the 
relevant <9-problem, and indicate how familiar steps for the Hermitian 
model, e.g. the g-function "undressing", might look like in the case of 
the normal model. We use the particular model considered recently by 
P. Elbau and G. Felder as a case study. 

1. Introduction 

In these notes we attempt to develop for the normal matrix model a 
formalism analogous to the Riemann-Hilbert method in the theory of Her- 
mitian matrix model. As in the latter case, the starting point is proper 
analytical characterization of the relevant orthogonal polynomials. Unlike 
the Hermitian matrix model, the orthogonality condition for the polynomials 
associated with the normal model is formulated with respect to a measure 
on the plane. This, as we will see below, leads to the replacement of the 
Riemann-Hilbert problem of [4] by a certain 9-problem. We shall present 
in detail the setting of the 9-problem for the case of what we will call in 
these notes the Elbau-Felder model. This model arises as a natural regu- 
larization of the normal matrix model of P. Wiegmann and A. Zabrodin 
in [H] by restricting the matrix integral of the latter to normal matrices 
whose eigenvalues lie in a compact domain D of the complex plane. Using 
the Elbau-Felder model as a case study, we shall also outline a possible in- 
version of the Deift-Kriecherbauer-McLaughlin-Venakides-Zhou (DKMVZ) 
asymptotic model. The DKMVZ method proved to be very efficient in the 
asymptotic analysis of the oscillatory Riemann-Hilbert problems appearing 
in Hermitian matrix model. We have not yet succeeded in providing com- 
plete generalization of the DKMVZ scheme for the orthogonal polynomials 
on the plane; in fact, we rather highlighted the challenging difficulties to 
be overcome. We hope, however, that these notes might stimulate further 
development of the analog DKMVZ asymptotic method for the orthogonal 
polynomials on the plane and related normal matrix models. 
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2. Preliminaries 

2.1. Normal matrix models and orthogonal polynomials. Let D be 

a bounded domain on the complex plane C containing the origin, and let 
V(z) be a real-valued smooth function on C. Following P. Elbau and G. 
Felder [3], we shall study the normal matrtix model characterized by the 
partition function Zn defined by the following TV-fold integral, 




Let xd be the characteristic function of the domain D. 

Definition 1. Orthogonal polynomials on C with respect to the measure 
e~ NV ^XD(z)d 2 z are polynomials P n (z) = z n + a n -\ n z n ~ l + • • • + ao n , sat- 
isfying 

(2.1) f I P n {z)P^(z)e- NV ^d 2 z = h n 5 mn for all m, n = 0, 1, 2 . . . . 

' D 

The following lemma is standard. 
Lemma 2. 

2V-1 

Z N = N\H h n . 

n=0 

The proof is exactly the same as in the case of the Hermitian model 
(see e.g. PQ). As in the case of the Hermitian model, Lemma [2] reduces 
the question of the asymptotic analysis of the partition function Zjy as 
N — s- oo to the asymptotic analysis of the orthogonal polynomials P n {z) as 
n,N — > oo. 



3. Matrix 5-problem 



Using the orthogonal polynomials on the line as an analogy (see [1], pQ), 
we set 



(3.1) 



YJz) 



I nE^l e -NV(z') d 2 z , 

7T J J Z'—Z 

D 



\ D 



It follows from the formula 



d 1 



- = 7t5(z - z'), 



dz z — z' 

understood in the distributional sense, that 

d 



(3.2) 



dz 



Y n (z) = Y n (z)(I-G(z)), 



ORTHOGONAL POLYNOMIALS AND d-PROBLEM 

where / is 2 x 2 identity matrix and 



(3.3) G(z) 



1 e- NV ^ X D(z) 
1 



The following proposition is central (cf. the case of the orthogonal polyno- 
mials on the line). 

Proposition 3. The matrix Y n {z) is the unique solution of the d-problem 
(|3.2p - (|3.3p with the normalization 



(3.4) Yn (z) = [I + 0[-))[ Q ^ n ) = (I + 



as \z\ — > oo, where 03 = (o -i)- 

Proof. It follows from the geometric series expansion 

oo / i\ k 



as \z\ — » oo, and the property (|2.ip . rewritten as 
(3.5) / / P n (z)z m e- NV ^d 2 z = h n 8 mn , 



D 



that the matrices (|3.ip satisfy normalization (|3.4p . 

Conversely, suppose that the matrix solves the 3-problem (|3.2p with 
the asymptotics (|3.4p . It follows from the special form f|3.3|) of the matrix 
G(z) that (Y n ) n (z) = P n (z) and (1^)21 (z) = Qn-i(z) — polynomials of 
orders n and n — 1 respectively, and 

J^«)l2(*) = -P^)e- NV ^ XD (z), 
^z(Y n ) 22 (z) = -Q^(z)e- NV ^x D (z). 

Now it follows from normalization (|3.4p that the leading coefficient of poly- 
nomial P„(z) is 1 and 



7T 



7T ./ ./ z' - Z 



D 
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Using geometric series and normalization (13.4)) once again, we obtain that 
polynomials P n and Q n -i satisfy 



P n {z)z m e- NV( - z) d 2 z = for m < n, 
Q n ^ 1 (z)z m - 1 e- NV ^d 2 z = 7r5 mn for m < n. 



D 



D 

From here it follows that 



P n {z)P m (z)e- NV{z) d 2 z = for all m < n, 



D 



which is sufficient to conclude that P n (z) are orthogonal polynomials on C 
with the weight e~ NV ^xo(z)- Finally, polynomials Q n (z) satisfy 

jj Qn-i(z)Pm-i(z)e~ NV( '^d 2 z = -ir5 mn for all m < n, 
D 

so that Q n {z) = ——P n {z). □ 

Similar to the case of the usual orthogonal polynomials, Proposition [3] 
reduces the asymptotic analysis of the orthogonal polynomials (|2.1|) to the 
asymptotic analysis of the solution of the 3-problem 



4. ELBAU-FELDER POTENTIAL. TOWARDS A NORMAL MATRIX VERSION 
OF THE DKMVZ ASYMPTOTIC APPROACH. 

We will consider the matrix model with the weight e~ NV ^xo{z), where 
V(z) is the Elbau-Felder [3] potential 

(4.1) V(z) = ^\z\ 2 -2Re^t k z k Y 

where t% = 0, \t%\ < 1/2 and toV(z) is positive on D\ {0}. Using again the 
Hermitian matrix model analogy, we shall expect that a fundamental role 
in the asymptotic analysis of the 3-problem (|3,2p - (|3,4p will be played by the 
equilibrium measure. 

Definition 4. An equilibrium measure for V on D is a Borel probability 
measure \i on D without point masses so that 

I(fi) = inf I(u), uCM(D), 

where M(-D) is the set of all Borel probability measures /i on fl without 
point masses, and the functional I{y) is defined by the equation, 

/(i/):= IV{z)dv{z)+ I I \og\z-C\~ l dv{z)dv{Q. 

J J Jz^C 
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Theorem 5 (Elbau-Felder [3j). There is 5 > such that for all < t < 5 

the unique equilibrium measure d[i exists and is given by 

d Kz) = -^-XD+(z)d 2 z, 
where the domain D + C D contains the origin and has the property that 



1 



*o = ~JJ d 2 z, 

D+ 

t k = / / z~ k d 2 z = i> zz~ k dz, k = 1, . . . , n + 1, 

irk J J 2irik J 

C\D + dD + 

t k = 0, j > n + 1. 

These relations determines D + uniquely. In fact, the boundary V of D + is 
a polynomial curve of degree n, i.e. V is a smooth simple closed curve in the 
complex plane with a parametrization h : S 1 C C — > C of the form, 

h(w) = rw + ao + a\w^ x + ... + a n w~ n , \w\ = 1, 

with r > and a n ^ 0. The equilibrium measure has the following properties. 
Set 

E{z) = V(z) + 2 fflog\z - Cr^MO- 



1. E(z) = Eq — a constant — for z € D + . 

2. E(z) >E — for z£ D\ D+. 

4.1. A naive DKMVZ scheme. Suppose that there is an analytic func- 
tion g(z) with the following properties. 

1. g(z) = log z + O (-) as \z\ — > oo. 

2. V(z) - g(z)- Wj_ = E on D+. 

3. V{z) - g(z) - g{z) > E on C\D + . 

Such function g(z) could be used to study the asymptotics of the matrix 
Y n (z) in the limit, 

n 

(4.2) n, N — > oo, — = 7, 7 is fixed, 

in exactly the same manner is it is done in the case of the orthogonal poly- 
nomials in the line (see [2] and [T]). Namely, set 

V 7 (z) = -V(z) 
7 

and consider corresponding equilibrium measure dfi^ (assuming that < 
7io < S) with the domain 1)4.(7) an d the corresponding function g^(z) 
satisfying properties 1—3 (with D + and Eq replaced by D + {^) and £"0(7)). 
Then we can "undress" the <9-problem (|3.2p - (|3.3p with the normalization 
flUJ) by setting 



nEn('y) f s. . nEnOy) 

— F^svn I ^\ ( ,ng- 1 {z)u 3 + — f±L 



Y n {z) = e —^ n (z)e 
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The resulting matirx *$> n (z) satisfies the simplified d-problem 



(4.3) = *„(^ ' 



with the standard normalization 

(4.4) ® n (z)=I + o(^\ as \z\ -> oo, 

which follows from property 1 of the function g~f{z). 

It is easy to pass to limit (14.2j) in the <9-problem (|4.3p ^ (|4.4p . Indeed, it 
follows from properties 2-3 of the function g-y(z) that 

lim V n (z) = ^o(z), 

n,N— >oo 

where the matrix ^>o(z) satisfies the following model d-problem. 

(4.5) JU„(*)^{P *™ 
with the standard normalization 

(4.6) V (z) = 1 + ( - J as |z|^oo. 



This model 5-problem is easily solved explicitly, 

, , A i si ^ 

*o(*) = D+(7) 

V° 1 

4.2. Function <?(.z). Of course, the main assumption that there is an an- 
alytic function function g(z) satisfying the properties 1-3 is not correct. 
Firstly, it follows from the property 1 that g(z) in the neighborhood of in- 
finity is defined up to an integer multiple of 2iri, which is not a drawback 
since e" 9 ^- )<T3 si well-defined for integer n. Secondly, the property 2 implies 
that the function V{z) is harmonic in D + , which clearly contradicts (|4.ip . 
Adding to the confusion is the formal manipulation 

log \z - C| 2 = k)g(z - C) + log(£ - 0: 

which suggests that 

(4.7) g(z) = jj \og{z - C)dM(C) = jfj j lo ^ z ~ 

D+ 

satisfies properties 1-3. However, this is not so since we need to treat 
carefully the branches of log in order to define the integral in (|4.7p and 
investigate its analytic properties. 
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For this aim, consider the logarithmic potential given by the uniform 
distribution of charges in the domain D, 

Vo(z) = J J log \z - w\ 2 d 2 w. 

D 

Let T = dD with fixed point Co £ T. For z € D denote by D e (z) domain 
obtained by removing the disk of radius e around z, so that 

dD £ (z) =TU-C £ (z), 

where C e (z) is the circle \w — z\ = e oriented counter-closkwise, and the 
minus sign denotes negative orientation. Since 

log \w — z\ 2 dw A dw = — dflog \w — z\ 2 wdw ) — = dw A dw, 

w — z 

by Stokes' theorem, we have 

Vo( z ) = ^ J J log \z — w\ 2 dw A dw 

D 

= - lim / / ( — d(los \w — z\ 2 wdw) dw A dw ) 

2e^0jJ V V ' ' ' W-Z ) 

D e (z) 

= — lim [id | log \w — z\ 2 wdw H — wdw 

2* e^O J J \ W-Z 

D e (z) 

dD e (z) 



Now 



lim / log |C - ^| 2 C^C = 0, lim / S-(d(=-2iri\z\ 2 
e^O J e^O J Q - z 



C e (z) C e (z) 

so that 

V (z) = tt\z\ 2 + 1 j> (log |C - z\ 2 QdQ + J— (d( ) . 
r 

Set lo = wdw and define a function Q on F \ {Co} by 

C 



n(c) = y u>, 



Co 

where the integration is along the oriented path in V connecting points Co 
and C- We have Sl_ (Co) = for a path consisting of a single point Co> and 

+ (Co) = j> wdw = jj dw Adw = 2iA(D), 
r D 
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for the loop T starting and ending at Co 5 where A(D) is the area of D. Thus 



j log |c - z\\dQ = j> log |c - z\ 2 dn(0 = A(io g \ z - Cl 2 o(0)£ 

r 



= 2iA(D) log \z-(o\ 2 -f n(C) + , 



so that 



V (z) = vr|z| 2 + A(2?) log \z-C \ 2 + \j dC + S§ <*C " ^ C<) ■ 

r 

Since the potential Vq is real-valued, we have 



V„(z) = 7r| 2 | z + /l(O) log \z - Co| z + \ f I ^ <i< - ^ <iC 



_ »(o-o(c) + ici 2 dC - \ 

Finally, observing that u; + u) = d|w| 2 , we get 

C 



n(C) + n(C) = / ^kl 2 = ICI 2 -|Col 2 , 



Co 

so that 



0(C) - 0(0 + |CI 2 = 20(C) + !Col 2 , 

and we obtain 

(4-8) 

V (z) = vr(|,| 2 - |Co| 2 ) + A(D) log \z-Q ? + \j ^ d C - ^ ^ , 

where z £ D. This is a desired representation of the area potential Vq(z) as 
the real part of the first derivative of a single layer potential. 
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The same computation for z € C \ D gives 



(4.9) V Q (z) = A(D) log \ z -C^+ l -j \W£L dC~^= dCj , 

Returning to Elbau-Felder potential V{z) and setting D = D+, T = dD + . 

we get E(z) = V(z) Vq(z), so that 

irt 

(4-10) 

v( z) - I (W - - icoi 2 ) - log i, - coi 2 - 4 / « - d c) = E 

when z E D + , and 



when z € D-. 

Remark 6. We note that equation (|4.10p . i.e. the statement that the l.h.s. 
of (|4.10p is constant when z € D + , is equivalent to the moment equations 
of Theorem [5] which determine the contour T (cf.[3], p. 12, Lemma 6.3). 



Now we are ready to introduce the function g(z). Namely, set 



(4.12) g(z) = log(* - Co) + ^— I ^ dC 

2vrt / C - z 



The function g(z) is holomorphic in C\T, is multi- valued with periods 27rzZ 
(singe- valued on the plane with the outside cut starting from Co) and has 
the asymptotics 

g{z) = log z + 0(z~ 1 ) as z — > oo. 

The function e n9<yZ ^ is single-valued for n € Z. The function <7(z) is discon- 
tinuous on r (by Sokhotski-Plemelj formula). 
We summarize this as the following statement. 

Proposition 7. The Elbau-Felder potential V(z) has the following repre- 
sentations 



(i) For z e D + , 

V(z 

(ii) For z G D_, 



V(z)-g{z)-g(z)=E + U\z\ 2 -U 2 ). 

to 



V(z) - g(z) - g{z) = E{z). 

(iii) For z G D\D+, 



(4.13) V(z) - g(z) - g(z) = E(z) > E Q . 
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4.3. A first possible version of the DKMVZ scheme. The correct 
strategy is now the following. Let g^{z), D+fa), £4,(7); etc. denote the 
respective objects associated with the potential Vy(z). We set 



ng-, (z)a 3 + — a 3 - CT3 



(4.14) y n (z) = e ^ CT3+ **> CT3 M/„(z)e' 
The resulting matirx \T/ n (z) satisfies the 9-problem (the correct version of 



3 T / \ 

oz 





v0 



-e 7 'o 




(4.15) 



— tt n (z) = * n (z) I -e 



n( E(z)-E (j)+ 





7*0 



XD+W , zGC\£> H 



= * n „(z)e^ Q{z)a3 , z g r = r( 7 ), 

with the standard normalization 



(4.16) 



*n(«) 



£+o( - 

Z 



as 2 



00. 



By virtue of condition (j4.13jl . we expect that the limiting function ^(z) 
satisfies the model problem 



(4.17) 



d_ 

0z 



-e t'o 




with the standard normalization 



^_(z)e^ Q(2)CT3 ! 



(4.18) 



* (z) = J + O ( - ) as \z\ 



z € U+( 7 ) 

zeT, 

— > 00. 



The open questions are now the following. 

(1) How to solve this model problem? The "unfortunate" thing is the 
presence of complex conjugation in (|4,17p . Indeed, if we neglect the 
jump across the contour T, the 9-problem alone can be of course 
solved explicitly, 



(4.19) 



*o(*) 



1 J // 
\0 1 



_n|C 

e "^b 72 



and the solution won't have any jumps. If not for the complex 
conjugation, the function ^o(z) could be used to undress in the 
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usual way problem ((HIT]) and reduce it to a pure Riemann-Hilbert 
problem. 

(2) The arguments that led us to the model problem (14. 17ft and which 
are based on inequality (14. 13f) . even on the formal level, are not 
very convincing: the real part of Q(z) is \z\ 2 — |£o| 2 ^ so that the 
diagonal jump matrix on T is not pure oscillatory. 

4.4. A second possible version of the DKMVZ scheme. The above 
deficiency of the proposed analog of the DKMVZ scheme can be partially 
overcome by performing the following modification. Let us replace the func- 
tion 0(0 by the function, 

1 (< 

Oo(C) = ~ {wdw — wdw). 

2 J Co 

The function Oo(C) is pure imaginary on T and it is related to the function 
0(C) by the equation, 

o (C) = «(C) - IC|2 ~ 2 ICo ' 2 - 

By using again Stokes' theorem, we observe that 

d( d( 



(ICI 2 -|Col 2 ) 



into J \C — z C — z 



1 ff ( C , C A ,2a i, uJ 



2vrt J Jd + VC - z Q-zJ *o 
This allows to re-write equations (|4.1U|) and (|4.11|) in the form, 



V(z)-lo g \z-Co\ -—fl—dC-j-jdt 



v 



(4-20) + -V // ( -i- + J-_ ) d 2 ( = E 



2vrt J Jd+ VC - z C, - z 

when z G D + , and 



when z € D_. These formulae in turn yield the following modification of 
the definition (|4.12|) of the (^-function 

• I fi oK) Jr i ff C ,2, 



(4.22) = log( , - Co) + _ « - _ yy ^ c . 



12 ALEXANDER R. ITS AND LEON A. TAKHTAJAN 

Note that the function g(z) is not holomorphic in C \ V anymore ! In fadO, 

(4.23) | 9(2) = _L XD+(Z) . 

A slight modification is also needed in the definition of the function ^ n {z)\ 
indeed, we should put, 

(4.24) Y n {z) =e 2 a *V n ( z )e n *'M ff * + —9r Lff3 . 

Taking into account (|4.23p . the 9-problem for the matrix now reads, 

(4.25) A^) + n _£_^) . 3= i^i)^ -A, z€D+ , 

(4.26) -V n (z) = y n (z)r Q Q X° + W\ zeC\D + , 

* n+ (z) = y n „(z)e^ n ° {z)as , z G r = T( 7 ), 

with the standard normalization 

(4.27) ty n (z) = I + O ( -) as \z\ -»■ oo. 



The function f2o(z) is now purely imaginary. Therefore, the arguments 
based on inequality (|4.13p seem to be more sound than in the previous 
approach, and they lead us to the following new model 9-problem 

(4.28) jL*o (z)+n ^o (z)CT3= ^^0 -1^ 

(4.29) JU#(*) = 0, zeC \D + , 

with the standard normalization 

(4.30) *J(2)=7 + oQ as |*| -»■ oo. 

Remark 8. Due to the presence of the large parameter n in the left hand 
side of equation (|4.25j) . the transition to the model problem (|4.28p - (|4.30p is 
still not quite satisfactory even on the formal level. 



Probably this property of the function g(z) reflects a major difference between the 
Riemann-Hilbert problem and the 9-problem. 
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4.5. An important concluding remark. In context of the theory of or- 
thogonal polynomials, a matrix 5-problem has also appeared in the recent 
work of K. McLaughlin and P. Miller [6] devoted to orthogonal polynomials 
on the unite circle with the non-analytic weights. However, unlike the prob- 
lem (I3.2p - (|3.4j) . the 5-problem of [6] is not the starting point of the analysis; 
indeed, the staring point of [6] is still the usual matrix Riemann-Hilbert 
problem and the <9-problem of McLaughlin and Miller is introduced out of 
the necessity to modify the "openning lenses" step of the usual DKMVZ 
scheme. Even more important difference between the 5-problem considered 
here and the 9-problem in [6] is the absence of the complex conjugation in 
the basic 9-relation. In one hand, this fact simplifies the implementation of 
the "undressing procedures" — the very important technical element of all 
integrable asymptotic schemes. On the other hand, as we have shown, the 
presence of the complex conjugation in the right hand side of (|3.2[) is truly 
essential for the incorporation into the asymptotic analysis of the concepts 
of equilibrium measure and (/-function. 

In spite of these differences we believe that using methods of [6j will help 
allow to overcome the indicated above obstacles in the asymptotic analysis 
of the <9-problem (|3.2p - (|3.4p . Specifically, we we think that one needs to 
develop and then apply to the model problems (|Q7|) - (|4TTg|) or (|4T28]) - (|4T30|) 
of the the 9-version of the "openning lenses" step in the DKMVZ method. 

In conclusion, we want to point out at the 9-method in the theory of 
integrable systems, introduced long ago by A.S. Fokas and M.J. Ablowitz 
in [5], as yet another source of tools for the analysis of the 9-problem (j3.2D - 

433}. 
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